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Abstract. Although fractional Brownian motion was not invented by 
Benoit Mandelbrot, it was he who recognized the importance of this 
random process and gave it the name by which it is known today. This 
is a personal account of the history behind fractional Brownian motion 
and some subsequent developments. 
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Since Benoit Mandelbrot's passing in October 
2010, many well-deserved tributes have been paid to 
him.^ Benoit influenced a great many fields ranging 
from the physical sciences to economics, and math- 
ematics was certainly among them. Benoit's great 
gift was his ability to recognize the hidden potential 
in certain mathematical objects.^ 

I had the good fortune to observe Benoit's math- 
ematical analysis in action, and I would like to tell 
you about my experience with one of the objects 
that Benoit worked with, the random process known 
as a fractional Brownian motion. Although fractional 
Brownian motion was introduced by Kolmogorov, 
it was Benoit Mandelbrot who recognized the rel- 
evance of this random process and, in his seminal 
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^There was a special symposium at the Ecole Polytech- 
nique in Paris in March 2011, one at Yale in April 2011 and 
a number of sessions related to Mandelbrot's work took place 
at the annual meeting of the American Mathematical Society 
in Boston in January 2012. 

^Benoit Mandelbrot studied with Paul Levy, who is widely 
acknowledged for his mastery of the Brownian world. 



paper with Van Ness [18], derived many important 
properties. There, he gave this process the name by 
which it is known today. See [15] for a general re- 
view. 

Let me recount first how I met Benoit. At the be- 
ginning of the seventies, I was a graduate student at 
Columbia University in the Department of Mathe- 
matical Statistics — a small department but home to 
prominent faculty such as Herbert Robbins, David 
Siegmund and Yuan Shih Chow. Although I had a 
fellowship during the academic year, I needed to find 
summer work — something I failed to do in my first 
year. I had sent my Curriculum Vitae to many com- 
panies in New York City, but I did not receive a 
single reply. 

For my second year, I decided to proceed differ- 
ently. I asked members of the Department for con- 
tacts. This is how I was put in touch with Benoit 
Mandelbrot, who was then at IBM Research — an 
hour's drive from New York City — but was also nom- 
inally an Adjunct Professor in the Department. In 
January of my second year, I called him and in- 
quired about potential summer jobs. The conversa- 
tion began in English but quickly turned to French. 
I had expected it to last a few minutes, but the con- 
versation lasted an hour with Benoit doing most of 
the talking (as was often the case). He ended the 
conversation, saying that he knew of no jobs. But 
a few months later, he called me back. As things 
developed, it turned out that he needed a program- 
mer for the summer and asked if I was interested. 
I accepted. This is how I became acquainted with 
his research at the time, which involved fractional 
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Brownian motion and its application to hydrology, 
and how I ended up as Mandelbrot's student. 

It started with the so-called "i^/S* statistic," where 
R is the range of partial sums of the data, and S 
is the sample standard deviation. It is a statistic 
that the British hydrologist Harold Edwin Hurst, 
in the first half of the twentieth century, had used 
to study the yearly variation of the levels of Nile 
river in Egypt [8]. The original work on the sub- 
ject by Benoit Mandelbrot appeared in 1965 in the 
Comptes Rendus [12]. 

Under the usual assumptions of finite variance 
and independent and identically distributed obser- 
vations, the R/S statistic should grow like n^/^, 
where n is the sample size. The Nile data, how- 
ever, indicated a growth of , where 1/2 < H < 1. 
The growth n^/^ is typically associated with ran- 
dom walk, so , with 1/2 < H < 1 must corre- 
spond to something else. This is why Mandelbrot 
suspected that a process like fractional Brownian 
motion Bfjit) may perhaps be relevant in this frame- 
work^ since, while the standard deviation of Brown- 
ian motion at time t is i^^^, that of fractional Brown- 
ian motion at time tist^, where < -ff < 1 [19]. The 
letter H, which refers to the hydrologist Hurst and 
which was used by Mandelbrot, has become stan- 
dard in this context, and it now labels the fractional 
Brownian motion. 

The term "fractional Brownian motion" was coined 
by Mandelbrot and Van Ness in the now classical pa- 
per [18]. Fractional Brownian motion has a number 
of nice properties, one of which is "self-similarity." 
A process {X(t),t G M} is self-similar with index 
H >0 if for any a > 0, the process {X{at),t G M} has 
the same finite-dimensional distributions as {a^ X(t), 
t G M}. Thus, like a fractal, there is scaling, but it 
is not the trajectories of the process that scale, but 
the probability distribution, the "odds." This is why 
this type of scaling is sometimes called "statistical 
self-similarity" or, more precisely, "statistical self- 
affinity." 

The fractional Brownian motion process is then 
characterized by the following three properties: 

(1) the process is Gaussian with zero mean; 

(2) it has stationary increments; 

(3) it is self-similar with index H, < H < 1. 

Fractional Brownian motion reduces to Brownian 
motion when H = 1/2, but in contrast to Brownian 



''Some hydrologists argue instead that the behavior may 
be due to nonstationarity. See [9] and [21]. 



motion, it has dependent increments when H ^ 1/2. 
Fractional Brownian motion was first introduced in 
1940 by Andrei Nikolaevich Kolmogorov [10], who 
was studying spiral curves in Hilbert space. It was 
considered by Richard Allen Hunt [7] in the context 
of random Fourier transforms and by Akiva Moisee- 
vich Yaglom [34] , who studied the correlation struc- 
ture of processes that have stationary nth order in- 
crements. However, it is undoubtedly the seminal 
paper of Mandelbrot and Van Ness which put the 
focus on fractional Brownian motion and gave it its 
name. Why the term "fractional?" This is because 
the process can be represented as an integral with 
respect to Brownian motion B{t), as follows: 

BHit)= f {{t-sr-'i'-{-sr-'i'}dB{s) 

(1) 

+ / {t-sf-^'^dB{s) 
Jo 

/oo 
{{t-s)'l-'/'-i-s)'l-'/'}dB{s). 
-oo 

The integrals are well defined because the integrands 
are square integrable with respect to Lebesgue mea- 
sure. The form of the integrands is also reminiscent 
of the one that appears in the n-fold iterated inte- 
gral formula, 

/ dtn-i / dtn-2- ■■ / dti / g{s)ds 
Jo Jo Jo Jo 

and therefore (1) can be regarded as involving "frac- 
tional integrals." This, in fact, turns out to be more 
than a superficial analogy! 

The focus on fractional Brownian motion has proved 
to be extremely fruitful because it has allowed all 
kind of extensions, some of which were hinted at by 
Benoit Mandelbrot. 

For example, the Gaussian noise "di?" in (1) can 
be replaced by an infinite variance Levy-stable noise, 
giving rise to the linear Levy fractional stable mo- 
tion, which is an infinite variance self-similar pro- 
cess with stationary, but dependent, increments [28]. 
The kernel can also be replaced by a random sum 
of pulses [3]. From a different perspective, the single 
integral in (1) can be replaced by a multiple inte- 
gral, so that it becomes an element of the so-called 
Wiener chaos [24, 30], of the form 

(3) / gt{xi,...,Xk)dB{xi)-- - dB{xk), 
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for a suitable kernel gt and where prime indicates 
that one does not integrate on the diagonals. More 
specifically, if one chooses 




then the resulting process (3) is also self-similar with 
index 1/2 < H <1 and has stationary increments. It 
reduces to fractional Brownian motion if A: = 1 but 
is non-Gaussian if > 2. The marginal distribution 
for = 2 is studied in [33] . 

The representation (3) with the kernel gt in (4) 
is called a "time representation," but there are also 
other representations, for example, a "spectral rep- 
resentation" or a "finite interval representation" [26] . 

One can also try to define stochastic integrals, 
where the integrator is dBn, even though fractional 
Brownian motion Bh does not have independent in- 
crements. One then needs to define integrals of the 
type 

(6) / g{x)dBH{x), 

first for nonrandom functions g [20, 25], and then 
for random functions g [2]. One can also consider 
stochastic differential equations driven by fractional 
Brownian motion [22]. 

In a more applied vein, one can focus on the in- 
crements 

(7) X{n)=BHin)-BHin-l), n > 1, 
which form a stationary time series with covariance 

(8) r{k) = E[X{0)X{k)] ~ Ck^^-'^ 

as /c — > oo. The Fourier transform of the covariance 
(spectral density), 

oo 

(9) /(A) = E '(^y 

blows up at the origin ii 1/2 < H < 1 since /(O) = 
X]^_o^r(0) = oo. This type of dependence is called 
long memory, long-range dependence or strong de- 
pendence [31]. Time series with long-memory are 
important in modeling, particularly in econometrics. 
Financial returns, for example, appear uncorrelated, 
but their squares often display long-memory [1, 4]. 

Mandelbrot was very interested in finance and in 
developing suitable models for financial returns [14]. 



Together with his students Adlai Fisher and Laurent 
Calvet [16], he introduced a multifractal model of as- 
sets returns, BniOit)), where Bh is fractional Brow- 
nian motion, and 9{t) is an independent multifractal 
process corresponding to "activity time." A multi- 
fractal process has stationary increments and satis- 
fies 

(10) E|^(t)|« = c(g)r('^)+\ t>0, 

where r(g) is not necessarily a linear function of 
q. Observe that fractional Brownian motion itself 
is multifractal, more precisely monofractal, because 
by self-similarity, E|5H(i)l^ = 

corresponding to 

the linear function r(g) = qH — 1. The idea of "sub- 
ordination," replacing "physical time" by "activity 
time" can already be found as early as 1967 in his pi- 
oneering paper with Howard Taylor, with Brownian 
motion instead of fractional Brownian motion [17]. 

Statistics also enters into the picture. Given a time 
series, how does one check that it displays long- 
memory? And if it does, how does one estimate HI 
There is quite a large literature on the subject [27]. 
Many of the available tests are graphical [32] or 
asymptotic in nature. The asymptotic ones are re- 
lated to central limit theorems but also to so-called 
noncentral limit theorems which arise as follows. 

Consider a time series {Xn,n G Z}, and let = 
Var(^^^-i^ Xfc). What is the limit of the normalized 
sum 

^ [nt] 

(11) — E^'^' 

k=l 

as n — > oo? It is typically Brownian motion in the 
case of weak dependence. But if n € Z} is long- 
range dependent, it could be Brownian motion, frac- 
tional Brownian motion or a non-Gaussian process. 
The study of limit theorems in this context is thus 
very important, and while there is a large litera- 
ture about it, there are still many open problems 
[5, 6, 11, 13, 23, 29, 30]. 

This story illustrates one of Benoit Mandelbrot's 
contributions to mathematics. He tended to focus on 
a concept or mathematical object whose importance 
was not recognized. He studied and developed it — 
at times rigorously, at times heuristically — with the 
frequent consequence that the object of his atten- 
tion became the basis of major subsequent develop- 
ments. Because he challenged accepted views, many 
of his ideas met with initial resistance. Ultimately 
though, Benoit Mandelbrot's influence on the course 
of mathematical thinking has been far-reaching. He 
will be greatly missed. 
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